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Abstract 

It is shown that information transmission inside a cell can occur by 
means of mechanical waves transmitted through DNA. The propaga- 
tion of the waves is strongly dependent on the shape of the DNA thus 
proteins that change the shape of DNA can alter signal transmission. 
The overall effect is a method of signal processing by DNA binding 
proteins that creates a "cellular communications network". 

The propagation of small amplitude disturbances through DNA 
is treated according to the mechanical theory of elastic rods. Ac- 
cording to the theory four types of mechanical waves affecting exten- 
sion (compression), twist, bend or shear can propagate through DNA. 
Each type of wave has unique characteristic properties. Disturbances 
affecting all wave types can propagate independently of each other. 
Using a linear approximation to the theory of motion of elastic rods, 
the dispersion of these waves is investigated. The phase velocities of 
the waves lies in the range 5 ~ 8A/ps using constants suitable for a 
description of DNA. The dispersion of all wave types of arbitrary wave 
length is investigated for a straight, twisted rod. 

Based on these findings, we propose all-atom numerical simulations 
of DNA to investigate the propagation of these waves as an alternative 
measure of the wave velocity and dispersion analysis. 
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Introduction 

Our treatment of DNA is based on the hypothesis that such a well or- 
ganized system as a cell must have a highly sophisticated system of com- 
munications. For instance, the cell cycle is the result of a series of complex 
self organized processes that must be well orchestrated. A model has been 
developed to demonstrate how these events may be initiated by critical con- 
centrations of specific proteins which shift equilibrium to favor the advance 
of the cell cycle 0. However, it is also known that the cell-cycle can be dis- 
rupted if conditions are not suitable. Similar arguments apply to virtually all 
cellular processes and in order to achieve the required checks and balances a 
method of communication is necessary. 

Is there the possibility that information is transmitted through DNA elec- 
tromagnetically instead of mechanically? In such case DNA will function as 
a transmission line which requires total internal reflection (TIR) of the ra- 
diation within the DNA. To achieve TIR, the wavelength of the radiation 
must be 5-10 times less than the diameter of the transmission line. Since the 
diameter of the DNA is close to 20 A the radiation must have a wavelength 
close to 2 A. This wavelength is close to atomic dimensions so diffraction 
should dominate. Furthermore, the energy associated with this wavelength 
radiation is on the order of 10 5 kcal/mol which is sufficient to destroy chem- 
ical bonds and therefore not easily managed biologically. For these reasons, 
we believe that the communication will be achieved mechanically. 

The remainder of the paper will demonstrate that the mechanical prop- 
erties of DNA have the necessary time and spatial dimension to support the 
propagation of information while interactions between proteins and DNA 
provide a mechanism by which this information is processed. 

This manuscript is based on a continuous medium model of DNA. Such 
an approach is well known and widely used to describe solids, liquids and 
gases (see, for example [0, |3|). Thus, when we speak of infinitely small ele- 
ments of volume, we shall always mean those which are "physically" infinitely 
small, i.e. very small compared with wavelength or radius of curvature under 
consideration, but large compared with the distance between the atoms in 
DNA. 
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1 System of equation for the elastic rod dy- 
namics 



We utilize an elastic rod model described by ||J] that has been parameterized 
to represent DNA ||. The equations of motion follows: 



/•l+fixf/..) = b . d[Cl - + a x (c ■ (n - Op 

at v y as v v 

+ fx (C-(f-fo)) (2) 

In this system of equations, equation (p represents the balance of force 
and linear momentum and equation (0) represents the balance of torque 
and angular momentum according to Newton's Laws. Equations (|3|) and (Q) 
are continuity relations expressed in a non-inertial reference frame. Here, 
s and t are independent variables representing time and fiducial arclength, 
respectively f\ The functions that we wish to study are the four three-vector 
functions Q, T, u, and 7. 

The matrices /, C and D and the scalar p are as follows: Matrix I is the 
linear density of the moment of inertia tensor. Matrices C and D represent 
the elastic properties of the rod (C3: Young's modulus, C\ and C2: shear 
modulus, -D3: torsional rigidity, D\ and D2: bend stiffness), p is the linear 
mass density and is equal to 3.22 x 10 _15 K/M for an isotropic model of DNA. 
The matrices have the following values for an isotropic model of DNA: 



1 Fiducial arc length as describe in m is the same as actual arc length only when the 
rod is neither extended nor compressed. During the propagation of a disturbance along 
the rod, the rod will generally experience extension in some places and compression in 
others so the actual arc length will generally be different from fiducial arc length. Actual 
arc length plays no role in the equations of motion 
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For isotropic bending C% = C% and shear D\ = D 2 we will denote C = C\ = 
C2 and D = Di = D 2 . These definitions will be used below. The above 
values have been adapted from |J, 0. 

— * 

In biological terminology the components of T correspond to the three 
DNA helical parameters describing translation T = (shift, slide, rise) and the 
components of Q correspond to the three DNA helical parameters describing 
rotation Q = (tilt, roll, twist). If we attach a local coordinate frame, denoted 
by {di,d2,ds}, to each base-pair the axes will point in the direction of the 
major groove, minor groove and along the DNA helical axis, respectively. The 
shape of the DNA can be described by a three-dimensional vector function 
r(s,t) which is related to T and Q by a suitable mathematical integration. 

In elastic rod terminology the three-dimensional vector function r(s, t) 
gives the centerline curve of the rod as a function of s and t, showing only 
how the rod bends. To show the twist, shear and extension of the rod we 
must attach "director" frames made of the orthogonal triples {d\,d2,ds} at 
regular intervals along the rod. The director frames are evenly spaced along 
the rod when the rod is not extended and are all parallel to each other 
(with g?3 pointing along the rod) when the rod is not bent or twisted. Any 
deformation of the rod will be indicated by a corresponding change in the 
orientation of the director frames. 

The vector 7 (r) is the translational velocity of r in time (space): 

dr ^ dr 
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dt ds 

The directors are of constant unit length, so the velocity of each director u 
(Cl) is always perpendicular to itself, and it must also be perpendicular to 
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the axis of rotation of the local frame. 



dd k _ -> dd h -* -> 
— = u x d k — — = llx d k 
at as 

We will denote V of the unstrained state by IV In general an unstrained 
rod can have any shape. The most simple case is when the rod is straight 
with unit extension, Tq = {0, 0, 1}, because in the unstrained state: 

§ = {0,0,!} 

Similarly, an unstressed elastic rod may have an intrinsic bend and/or 
twist, denoted by flo, which must be subtracted from the total fl to give the 
dynamically active part. The simplest case is a rod with no intrinsic bend 
or twist Q = {0, 0, 0} and f = {0, 0, 1}. 

For short wave lengths (Section 0) we will suppose that fl = 
(floi, flo2, ^03) = const and T = (IV, IV, IV ) = const. In section |3] 
we will consider a straight rod with an intrinsic twist Q = {0, 0, fio} and 

f = {o,o,r }. 



2 System for small amplitude waves 

The results from this section are valid for any shape in which the curvature 
of the rod is much greater than the wavelength of the disturbance being 
propagated regardless of the intrinsic shape of the rod. 

Let us search as usual for the equilibrium point of the system ([[]-§) by 
setting 7, T, uj and fl to constants. It is easy to see that an equilibrium point 
of equations ([[]-§) is7 = cJ = 0, r = r and fl = Qo 

We shall suppose that each variable differs slightly from the equilibrium 
value so we retain only the linear terms in the equations. The linearized 
system (|l]-[|) is: 



— + ^xr = ^ + fi„X7 (10) 

— * 

dQ ^ =t duj . . 

In the system (|8[-|TT|) and further in this paper 7, T, uj and Q denote the small 
deviations from the equilibrium values, e.g. T — To — * T. 

In the usual manner we search for a solution to (|8|-|ll|) using harmonic 
analysis. We assume that each variable depends on time and space as follows: 

G l (s,t) = G t0 -e-^ t + iks (12) 

where Gi(s,t) denotes one of the four vector variables from the system of 
equations (|8|-|TTD and G i0 = const denotes its amplitude.^ 
Let us substitute (fT^) to the system (|8T-|TT|): 



- iujp^ = ik(C ■ T) + Vl x [C ■ rj (13) 

-uj) = ik (D ■ Q\ + tt x (b ■ VL) + f x (d ■ f) (14) 

— iuf + wxf = ik^f + Q x 7 (15) 
— iu>Q + uj x Q = ikul (16) 

For short wavelength we search for solutions with u oc k, so in the limit 
■ — > 00 we neglect all terms which contain cross products and the system 
ll3|-|i~6|) becomes: 



{uo 2 p-k 2 C)f ± = (17) 

(u 2 p-k 2 C 3 )T 3 = (18) 

{u 2 I-k 2 D)Q L = (19) 

{2u 2 I - k 2 D 3 )n 3 = (20) 

Here division into longitudinal and transversal components has been obtained 
by defining r_|_ = (r 1; T 2 , 0) and Q± = (0 1; Q 2 , 0). 

Equations (p!7|-^0|) represent four types of waves that can propagate in the 
elastic rod with velocity, V (according to the order of equations above): 

2 Note that a; is a scalar quantity corresponding to a frequency of oscillation and that 
to is a vector quantity that describes a rotation, so there should be no confusion between 
the two variables. 



6 



Shear waves (TjJ: 



8.16 ■ lO" 10 [KM/s 2 
1 ■ s ' / " '"• s ' 7 ~ \ 3.22-10- 15 [K/M] 

Extension waves ^3): 

^Extension 



5.03A/ps (21) 
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p \ 3.22- 10- 15 [K/M] 



8.2A/ps (22) 



Bend waves 



V, 



Bend 



2.7- 10~ 2S [KMys 2 
V 4.03 ■ 10- 34 [KM] 



8.2A/ps 



(23) 



Twist waves (fta 



Twist 




2.06- 10- 28 [fTM 3 /s 2 ] 
V 2-4.03- 10- 34 [iTM] 



5.0A/ps (24) 



These results were obtained for a rod with arbitrary shape because all 
terms that define the rod shape in equations ( fn^|16"l) were omitted. So, if 
wavelength tends to zero| (is the least space parameter of the problem) four 
wave types can propagate along the rod. 

It is well known that the measurement of wave velocity is a usual method 
for determining elastic properties of solids. So, wave velocity (transversal and 
longitudinal) and Young's modulus (shear modulus) are uniquely determined 
by formulae similar to the ( |2"T| - [2"4] ) . 

Hakim et al. || have determined the velocity of sound in DNA. Based 
on elastic constants obtained from |J , and used in this work, the velocity 
is approximately two times less than that obtained by Hakim et al. p. 

Experiments to measure elastic properties of single DNA molecules have 
been reported using scanning tunnelling microscopy 0, fluorescence mi- 
croscopy [0, fluorescence correlation spectroscopy [O], optical tweezers 



3 We must note once more that "wavelength tends to zero" means that it is still much 
greater than distance between atoms, see Introduction). 
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r^][jTB|, bead techniques in magnetic fields [14j JT5[, optical microfibers [|T| 
low energy electron point sources (electron holography) [|l7j] and atomic force 
microscopy (AFM) |T8|] . Each method differs in the molecular properties 
probed, spatial and temporal resolution, molecular sensitivity and working 
environment, so each method gives close but not the same value of elastic 
constants. So if we evaluate extension/compression velocity using constants 
obtained from the data reported in [12|] the result will be two times slower 
than the value reported here and based on data from and [[?]. 



3 Linear waves in the straight rod with no 
intrinsic bend or shear. 

Here we concentrate our attention on the case of the straight, twisted rod with 
no intrinsic bend or shear. In this section we suppose that Qq — (0, 0, Qq) 
and r = (0,0, To). As we will see this case allows for a complete analyti- 
cal analysis of extension/compression (Section jOl) , twist (Section |3~2] ), and 
bend/shear (Section j3.3|) waves of arbitrary wave length. Equations (||-|ll|) 
in component form are: 
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(32) 
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(33) 
(34) 
(35) 
(36) 

It is easy to see that equations (p7|) and (|33| ) (also ( |30| ) and (|36|) ) are 
independent from all other equations and describe extension (twist) waves. 
These two wave types will be discussed in the next two sections. 

3.1 Sound (extension/compression) waves 

From equations (|27]) and (|33|) it is easy to obtain two wave equations for the 
extension/compression waves: 

- < 37 » 

# " ^ = (38) 

These small amplitude waves have velocity ( p2|) and these equations have an 
harmonic solutionF]: 

4 Equations (]37[-p8[) have not only harmonic solutions. In general the solutions will be 
(d'Alembertian waves): 

73=Pi(s-*Wy)+P2(s + iy^) (39) 
r 3 = qi ( S -J^-) + q 2 (s + J^-) (40) 



where pi 7 2 and are arbitrary functions andpi,2 = =F\/C 3 /p <7i,2- The point to keep in 
mind is that the amplitude of the sound waves must be small. It is easy to see the movement 
follows directly from equations (^9|) or ( f40|). It i s clear that the arbitrary function pi has 
the same value for all arguments s — t^JCjJp = const. This relationship describes a 
straight line in the (s, £)-plane so pi is constant along this line. The same statement holds 
for all neighboring points so an initial shape moves along the s — t^JC^/p = const line with 
velocity \fCzJp- The arbitrary function p 2 describes motion in the opposite direction. 
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/ ^3 / ^~^3 

T 3 = A\ sin(s — tJ — ) + A 2 cos(s — tJ — ) 

73 = Bi sin(s — t\ — ) + B 2 cos(s — tJ — ) 
VP VP 



(41) 
(42) 



where A\^ are arbitrary constants (wave amplitudes) and B\ y2 = —\jC^/p A\^- 

Using elastic constants obtained from data in |J and the velocity of 
sound in DNA is equal to V soun d = 8.2A/ps (see equation fl22|)) and dis- 
persion is linear. This is the velocity of propagation of harmonic waves in 
DNA according to the linear approximation. In the linear approximation the 
amplitude of the wave must be small, but it can have any wavelength. 



3.2 Twist waves 

In the same manner equations 
describe twist waves: 



and (|36|) yield two wave equations which 
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21- 



dt 2 

d 2 n, 



- D 



D 



ds 2 











dt 2 ds 2 
Linear solutions of these equations can be written as: 

+ P 2 cos(s 





—) + Q 2 cos(s -tW— ) 

where P\p are the arbitrary constants (wave amplitudes) and Q\ 2 
tD 3 /2I P 1>2 . 



(43) 
(44) 

(45) 
(46) 



These solutions describe twist waves with velocity V tW i S t = yD 3 /2I rs 
5.0A/ps (see equation (p4|)). In the general case, the solutions of these equa- 



tions will be d'Alembertian waves described by (p9|-^0|). Again the dispersion 
law is linear. 
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3.3 Bend and Shear Waves 

Substituting ( |T2"D into the remaining equations of system ( p5[]36| ) yield: 



—iuICu± 
-iuT± + ujj_ x T 



— iuQi + 



ikcr ± + c n x r ± 

ikDQ ± + D p x Q_ 
ik^± + fi x 7 ± 
iku)j_ 



+ C 



TnXT 



(47) 
(48) 
(49) 
(50) 



Here, the definitions 7j_ = (71,72,0), Tj_ = (ri,r 2 ,0), uj± = ((Vi,co 2 ,0) 
and Qj_ = f^2 ? 0) are used. This is an homogeneous system of linear 
equations with unknowns 71, and Q±. Solutions other than the 

trivial solution exist only if the determinant of the coefficients is zero. This 
condition is satisfied if: 



{u 2 i - (k 2 + n 2 )Dy - 4k 2 D 2 ni 



2kpDcu 2 CT (Qq - r, 



2kCVL + - 

{u 2 I -{k 2 + Q 2 )Df -Ak 2 D 2 Q 4 0i 

Let us discuss the untwisted rod (Qq = 0)0. Equation (|5l"D yields: 



(51) 



u, 2 p - k 2 C - 



P lu 2 CT 2 







(52) 



oj 2 I - k 2 D 

It will be convenient to use dimensionless variables (marked by asterisk) 
uj* = u\Jl/C, k* = k^Jl/p in this section and define G 2 = ^!j- p . Equation 
( |52|) becomes(asterisks are omitted): 



(uo 2 -k 2 )(uj 2 -G 2 k 2 ) 



V 2 , , 2 



with solutions: 

^1,2,3,4 



±lyiG + l} 2 k 2 + T 2 ±^[G-l} 2 k 2 + T 2 } 



(53) 



(54) 



5 For untwisted rod system (J23-^6|) consists of two independent subsystems. Equations 
(J25|, p9|. p2, p5|) and (Eq, E8l p2, p4|) present two different polarizations of bend/shear wave. 
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Figure 1: Dispersion of Bend/Shear Waves in DNA. 




There are a total of four solutions. The ± sign before the braces merely 
changes the direction of propagation and will not be discussed further. The 
± sign between the roots produces two different branches of solutions as 
presented in Fig.[L| The "upper" branchP] corresponds to bend waves and the 
"lower" branch to shear waves. 

For a rod of infinite length, any value of u can be chosen that satis- 
fies equation fl53|), but for a rod of finite length there will be an additional 
constraint that must be satisfied. If the rod is fixed at both ends then the 
usual restriction of having a node at each end applies. This introduces the 
constraint L = (n/2)X between the length of the rod, L, and the possible 
wavelengths, A, so that only discrete points along either sub-branch in Figure 
1 will be observed. 



4 Conclusion. 

The primary results are that, in the linear approximation, four different types 
of waves can propagate through a uniform elastic rod. These waves corre- 

In this way we will name branch that does not go to zero in the limit k — > 0. 
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spond to extension(compression), twist, bend and shear. An extension or 
twist wave will propagate without exciting other modes or changing shape, 
and it has a linear dispersion relation. Bend and shear waves behave rather 
differently and have a nonlinear dispersion relation. Each type obeys a dis- 
persion law that describes two additional sub-branches. 

Utilizing constants suitable for DNA we find that, in the limit of small 
wavelengths, extension and bend propagate with a velocity of approximately 
8A/ps and twist and shear propagate with a velocity of approximately 5A/ps. 

Is is also significant that the dispersion relation for bend and shear waves 
is coupled with the inherent twist of DNA. We propose that these physi- 
cal phenomena enable proteins interacting with DNA to accomplish highly 
sophisticated tasks. For example the difference in extension and twist ve- 
locities can be utilized to measure the distance between two points on the 
DNA. Since the dispersion law for bend and shear waves depends on intrin- 
sic twist, a mechanism for measuring DNA topology exists because of the 
relation between twist, linking number and writhe. Other protein-protein 
communications can certainly be established to assist cellular mechanisms. 

4.1 Suggestion for a simulation experiment 

We suggest a molecular dynamics simulation experiment to check the corre- 
spondence between the linear theory of the propagation of waves in elastic 
rods and an all atom simulation of DNA. In all-atom molecular dynamics 
simulations of DNA the number of base pairs that can be simulated for a 
significant length of time is from tens to less than hundreds of base pairs. 
For a simulation of DNA with fixed ends it is sufficient to apply a sharp im- 
pulse (5 function excitation) to one end of the DNA and measure the time of 
propagation of this disturbance to the other end as a measure of the velocity 
of extension/compression waves in DNA according to equation (|22|) . The 
propagation of twist can be similarly tested by applying a torque. 

To measure dispersion one must make some simple calculations before the 
simulation and then excite one end of the DNA with a driving force of the 
appropriate frequency. In this manner a standing wave can be established in 
the DNA. The wave frequency and wave number (2n/X) are related by the 
well known relation: 

^ = V Wave k (55) 
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where u = 2irf is the wave frequency, Vwave is the wave velocity and k = 
2n/\ is the wavelength. In this case the wave length (wave number) is easy 
to evaluate: 

L = n—, where n— 1,2,3,.. (56) 

and / is the frequency of the driving oscillation. The values obtained from 
simulation can be compared with the number obtained from formula (^ or 
|53| ) as a means of correcting our choice of constants. 
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